Multiparameter entangled state engineering using adaptive optics 
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We investigate how quantum coincidence interferometry is affected by a controllable manipulation 
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I. INTRODUCTION 

Quantum entanglement |l| is a valuable resource in 
many areas of quantum optics and quantum information 
processing. One of the most widespread techniques for 
generating entangled optical states is spontaneous para- 
metric downconversion (SPDC) 0, i, i, 0. SPDC is a 
second-order nonlinear optical process in which a pump 
photon is split into a pair of new photons with conserva- 
tion of energy and momentum. The phase-matching re- 
lation establishes conditions to have efficient energy con- 
version between the pump and the downconverted waves, 
called signal and idler. This condition sets also a specific 
relation between the frequency and the emission angle 
of down converted radiation. In other words, the quan- 
tum state emitted in the SPDC process cannot be fac- 
torized into separate frequency and wavevector compo- 
nents. This leads to several interesting effects where the 
manipulation of a spatial variable affects the shape of the 
polarization-temporal interference pattern. For example, 
the dependence of polarization-temporal interference on 
the selection of collected wavevectors was studied in de- 
tail in [6!]. 

Here we engineer the quantum state in the space of 
transverse momentum and we study how this spatial 
modulation is transferred to the polarization-spectral do- 
main by means of quantum interferometry. We will fo- 
cus on type-II SPDC using birefringent phase-matching 
since the correlations between wave- vectors and spectrum 
are stronger than employing other phase-matching con- 
ditions. 

Our aim is twofold. From one point of view, we study 
the effect of spatial modulations on temporal quantum 
interference. This could be useful, for example, in quan- 
tum optical coherence tomography (QOCT) PjISI- When 
focusing light on a sample with non-planar surface, the 
photons will acquire a spatial phase distribution in the 
far-field, which may perturb the shape of the interference 
dip. Our results will provide a tool to understand this 
effect. 

From a second point of view, we would like to study 



and characterize spatial modulation as a tool for quan- 
tum state engineering. This may find application in the 
field of quantum information processing, where it is im- 
portant to gain a high degree of control over the pro- 
duction of quantum entangled states entangled in one or 
more degrees of freedom (hyper-entanglement). 

We start (Section II) introducing a theoretical model of 
a type-II quantum interferometer, comprising the polar- 
ization, spectral and spatial degrees of freedom. A modu- 
lation in the wave- vector space is provided by an adaptive 
optical setup and equations for the polarization-temporal 
interference pattern in the coincidence rate are derived. 
In our analysis (Section III) we will highlight and dis- 
cuss theoretically two interesting special cases. The first 
one is the possibility of restoring high visibility in type- 
II quantum interference with large collection apertures. 
In some situations, to collect a higher photon flux or a 
broader photon bandwidth, it can be useful to enlarge 
the collection apertures of the optical system. But, when 
dealing with type-II SPDC in birefringent crystals, for 
large collection apertures the effect of spatial walkoff in- 
troduces distinguishability between the photons, leading 
to a reduced visibility of temporal and polarization quan- 
tum interference. We will show that high visibility can 
be restored with a linear phase shift along the vertical 
axis. 

The second effect is the spatial counterpart of spec- 
tral dispersion cancellation [l^] . In the limit of large 
detection apertures, the correlations between the pho- 
tons momenta will cancel out the effects of even-order 
aberrations, exactly as in the limit of slow detectors the 
frequency correlations cancel out the even-order terms of 
spectral dispersion. The experimental demonstration of 
this effect has been reported recently ll]. 

Finally, in Section IV, we introduce a numerical ap- 
proach for practical evaluation of the results of the the- 
oretical model, discussing a few examples. By means of 
this approximated model, we will examine under what 
conditions the even-order aberration cancellation effect 
can be observed (Section V). 
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II. THEORETICAL MODEL 



Therefore the electric field at the position r and time t 
can be written as: 
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FIG. 1: Scheme of the proposed setup. Horizontally-polarized 
photons from type-II SPDC are assigned a phase dependent on 
the photon transverse momentum (j)o{qo), while the vertically- 
polarized ones are assigned a phase (j>e{qe)- The modulated 
photons enter a type-II quantum interferometer, which records 
the coincidence count rate as a function of the delay t between 
the photons given by an appropriate delay-line. 

Consider the scheme in Fig. [1] A laser beam pumps a 
X^^^ nonlinear material phase-matched for type-II para- 
metric down-conversion, creating a pair of entangled pho- 
tons. Each of the generated photons passes through a 
Fourier transform system, then enters a modulation sys- 
tem which transforms each transverse wave- vector for the 
horizontally polarized photon according to the transfer 
function Gi(qi) and for the vertically-polarized photon 
according to G2(q2). After being modulated in the q- 
space, the photons enter a type-II interferometer. A 
non-polarizing beamsplitter creates polarization entan- 
glement from the polarization-correlated pair emitted by 
the source. The beams at the output ports of the beam- 
splitter are directed towards two single photon detectors. 
Two polarizers at 45 degrees before the BS restore in- 
distinguishability in the polarization degree of freedom. 
An adjustable delay-line r is scanned and the coincidence 
rate R{t) between the detection events of the two detec- 
tors is recorded. An aperture is placed before the beam- 
splitter to select an appropriate collection angle. 



A. Notation 

Consider a monochromatic plane wave of complex am- 
plitude E{r) = Eqc"'^^''' with r = {x,y,z). For a given 
wavelength A, corresponding to a frequency Q, the wave- 
vector can be split in a transverse component q = {qx, Qy) 
and a longitudinal component /3(q; 



k=[q,/3(q;r!)] 



The wave-number is: 



k{n) 



(1) 



(2) 



The longitudinal component of the wave- vector is: 

/3(q,l)) = Vfc2(f^)_|q|2 (3) 



where p — {x,y). 

In paraxial approximation |qp <C so that: 

(3iq,n) ^ kin) - 



(4) 



2k{n) 



(5) 



For a quasi- monochromatic wave-packet centered around 
the frequency Qq one can write Q ~ flo + v, with v <ti D,q 
and this expression can be approximated by: 



iqr 

uq 2ko 



(6) 



where fco = fc(rio) and uq = (^^^n^|n=Oo) is the group 
velocity for the propagation of the wave-packet through 
the material. 



B. State Generation 

Using first-order time-dependent perturbation theory 
the two-photon state at the output of the nonlinear cystal 
can be calculated as: 



W ~ -j: I dtHiit) |0) 



where the interaction Hamiltonian is: 



(7) 



Hi{t) = ^ j drx^^Hr)El+^{r,t)E(-Hr,t)Et\r,t) 

(8) 

The strong, undepleted pump beam can be treated clas- 
sically. Assuming a monochromatic plane-wave propa- 
gating along the z direction: 



Ep{v,t) = Epe 



i{kpZ—ujpt) 



(9) 



The signal and idler photons are described by the follow- 
ing quantum field operators: 

E[-\v,t)= j dq, j dc^,e^[^(i"'^-)-+i'-''-'-'*la(q,,c^,) 

(10) 

The biphoton quantum state at the output plane of 
the nonlinear crystal is [12l |: 



\-<P) ^ j dqj dv^q, v)dl{q, + iy)dl{-q, Qq - u) |0> 

(11) 

Two photons are emitted from the nonlinear crystal, one 
horizontally-polarized (ordinary photon) and the other 
vertically-polarized (extraordinary photon), with anti- 
correlated frequencies and emission directions. 
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In the case of a single bulk crystal of thickness L and 
constant nonlinearity Xoj the probability amplitude for 
having the signal photon in the mode (q, f^o + 1^) and the 
idler in the mode (— q, flo — v): 



$(q, ly) = sine 



LA(q,z.) 



(12) 



For type-II coUinear degenerate phase-matching, the 
phase- mismatch function A(q, v) can be approximated 
to be: 



A(q, ly) = -vD + A/ea • q 



2|qp 



(13) 



where D is the difference between the inverse of the group 
velocities of the ordinary and extraordinary photons in- 
side the birefringent crystal and the quadratic term in 
q is due to diffraction in paraxial approximation. The 
last term is the first-order approximation for the spatial 
walk-off. 



This probability amplitude represents the superposi- 
tion of two possible events leading to a coincidence count 
in the detectors: 

1. the V-polarized photon with momentum qe and fre- 
quency LJe going through the lower branch to arrive 
at position xa in detector A, while the H-polarized 
photon with momentum qo and frequency lOq goes 
through the upper branch to arrive at position xb 
in detector B. 

2. the V-polarized photon with momentum qe and fre- 
quency LOe going through the lower branch to arrive 
at position in detector B, while the H-polarized 
photon with momentum qo and frequency lOq goes 
through the upper branch to arrive at position xa 
in detector A. 

Since the superposition is coherent, there are quantum 
interference effects between the two probabilities ampli- 
tudes. 



C. Propagation 



State engineering section 



Consider a photon described by the operator Sj (q, fi) 
(polarization j — e, o, frequency fi, and transverse mo- 
mentum q). Its propagation through an optical system 
to a point x^. on the output plane is described by the 
optical transfer function ifj (x^ , q; 17) . In our setup, the 
field at the detector will be a superposition of contribu- 
tions from the ordinary and extraordinary photons. The 
quantized electric fields at the detector planes are: 

e'^+\-kaM) = JdqJ dtJe'-*^[iJe(xA,q;w)ae(q,w) 

-I- i/o(xA, q; t^)ao(q, w)] 

E^+\KB,tB)^ JdqJ dLje'-*«[i/e(xB,q;^)ae(q,w)+ 

i7o(xB,q;w)ao(q,w)] 

(14) 

The probability amplitude to detect a photon pair at 
the detector planes, with space-time coordinates (xa, Ia) 
and (xb,<b), is: 

A{^A,^Ti;tA,tB) = {0\E'^/\xA,tA)E'^B\^B,tB)\^) 

(15) 

For the biphoton probability amplitude we get: 



A{xA,XB;tA,tB) = J dqodq_edujoduje^{qo,<ie;^o,i^e 
ffe(xA,qe;^.)i?o(xB,qo;Wo)e^*('"=*-^+"°*^) + 
i/o(xA, qo; Cc;„)ffe(xB, qe; c^,)e-^('^°*^+"=*«)l 



In the state engineering section, each of the two 
branches consists of a pair of achromatic Fourier- 
transform systems coupled by a spatial light modulator 
or a deformable mirror. Each Fourier-transform system 
consists of a single lens of focal length /, separated from 
the optical elements before and after it by a distance /. 
The first Fourier system maps each incident transverse 
wave- vector q on the plane II/tv to a point x(q) on the 
Fourier plane 11^?: 



x(q) = T^q 



fen = 



(17) 



where f is the focal length of the Fourier-transform sys- 
tem. Since we assume the system is achromatic for a 
certain bandwidth around a central frequency ilo, the 
position x(q) depends only on q and not on w. 

The spatial modulator assigns a different amplitude 
and phase to the light incident on each point, as described 
by the function G(x) = t{:x.)eJ'^'^^\ Each point is then 
mapped back to a wave-vector on the plane Hqut by 
the second achromatic Fourier-transform system. 

Using the formalism of Fourier optics [3], the trans- 
fer function between the planes Hin and Hqut can be 
calculated to be: 



/ii(xi,X3) =y"dx G(x) e~'7'"("i+"^) 

The corresponding momentum transfer function is: 

/ 



(16) 



ifi(qi,q3) = G 



fco 



'5(qi - qa) 



(18) 



(19) 
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2. Interferometer 

After the plane ^out the two photons enter a type-II 
quantum interferometer. Each propagates in free space 
to a birefringent delay-hne and a detection aperture p(x) 
to be finaUy focused to the detection planes by means of 
lenses of focal length /p. Following the derivation in [Gj 
the transfer function is: 



ff2(xi,q;w) = y" /i(xi,Xi;w)e''i-^Mxi 



— g»('-j/c)(di+<i2+/o)gxp 



Xi 



2 /rf2 



-i(crfi/2a;)|qpp ' ^ 



2c/o V/o 
Xi - q 



^ - 1 



c/o 



(20) 



where P(q) is the Fourier transform of |p(x)p. 

A combination of the two different stages is described 
by the transfer function: 



coherence of the photons and their area is larger than the 
spot into which the photons are focused by the collection 
lens, we integrate over the spatial and temporal coordi- 
nates. Therefore the coincidence count-rate expressed in 
terms of the biphoton probability amplitude is: 



R{t) = J dXA J dxB J dtA J dtB\A{xA,XB;tA,tB)\'^ 

(22) 

Following the derivation described in Appendix A, one 
gets: 



R{t) = i?o - A ( 1 - ) Wair) (23) 



where A(x) is the triangular function: 



Ha (xj , q« ; tJct) — Go 



/ 



i72(xj,q„;c^), (21) 



where the two functions Gi(q) and G2(q) are the mo- 
mentum transfer function which describe the modulation 
imparted respectively on the ordinary and the extraordi- 
nary photon. 



D. Detection 

Since the single-photon detectors used in quantum op- 
tics experiments are slow with respect to the temporal 



A(x) 



1 - 
0, 



\x\<l 
\x\ > 1 



(24) 



Therefore, the coincidence count rate R{t) is given by 
the summation of a background level Ra and an inter- 
ference pattern given by the triangular dip, A (l — -g^), 
that one gets when working with narrow apertures, mod- 
ulated by the function Wm{t) which depends on the de- 
tails of the adaptive optical system. 

The expressions for Rq and Wm{t) are: 



J 



Ro = J dq J dq'smc[MLe2 ■ (q- q')]Gt (^£q^ Gi (^£q'^ G^ 



,Mi62-(q-q') 



-[|q|^-|qT 



fco 



G2 { --pq' 



PA(q-q')^^s(-q + q') 



(25) 



and 



WMir) 



dq J dq'sinc 



MLe2 ■ (q + q')A 1 - 



2t 



g; 



ko 



fco 



G2 ( -^q' 1 e"' 



fco' 



) Gi 1 -q 



fco 



i^re..(q-q')g»^[|q|^-|qT]p^ [q + q/] [_ ^ q.)] 
I 



(26) 



In the following we will assume there is spatial mod- 
ulation only on one of the photons, therefore having 
G2(q) = l. 



III. PARTICULAR CASES 

Let's examine Eq. [23] in a few simple cases. 
First we will consider the case when no spatial mod- 
ulation is assigned to the photons and Eq. [53] will re- 
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duce to the results already described in the literature 
for quantum interferometry with multiparametric entan- 
gled states from type-II downconversion 0]. Then we 
will examine the effect of a linear phase, describing its 
implications for the compensation of the spatial walk- 
off between the two photons. Finally we will describe 
what happens in the approximation of sufficiently large 
detection apertures, introducing the effect of even-order 
aberration cancellation. 



A. No phase modulation 

Applying no phase modulation, our equations reduce 
to the ones derived in [6|]. Particularly we find: 



and 
Wg{t) = Sine 



i?0 =PA(0)Pi3(0) 



A'PLkp / 2t 
2diD ^ \}^Ul 



(27) 



Pa 



Mk 



2d,D 



re2 



Pb 



Mk 
2d7D 



p 

Tea 



(28) 



The shape of the interference pattern is essentially given 
by the product of the triangular function by the Fourier- 
transform of the aperture function, centered at t = 0. 
For physically relevant parameters the sine function is 
almost flat in the region where the triangular function is 
not zero. 




100 



x(fs) 



200 



300 



FIG. 2: On the right side we can see the interference patterns 
with three different detector aperture sizes: the corresponding 
aperture functions are shown on the left side. 



In this case the solution is quite simple: 



R{t) = i?o 



1 - A 1 



DL J 



with: 



2diD 
kpMRc 



(30) 



(31) 



Typically, sharp circular apertures are used in exper- 
iments. In this case, the function P[q] is described in 
terms of the Bessel function Ji{x). For a circular aper- 
ture of radius R: 



Ji(2i?|q|) 
i?|q| 



(32) 



However the Gaussian approximation is still a good one 
if the width Rg of the Gaussian is taken to roughly fit 
the Bessel function (of width R): in our case we take 
Rg = R/ {2^/2). 

Therefore Eq. ([501) still approximately valid in the case 
of sharp circular apertures, just taking: 



A^/2dlD 
kpAIRs 



(33) 



Mathematically, in Eq. ((30)l the interference pattern is 
given by the multiplication of a triangular function cen- 
tered at T = DL/2 by a Gaussian function centered at 
T = 0. The width of the Gaussian function ti decreases 
with increasing radius of the aperture Rb- Therefore, 
in the small aperture approximation, the width of the 
Gaussian is so large that it is approximately constant 
between t ~ and r = DL/2, giving the typical trian- 
gular dip found in quantum interference experiments. On 
the other hand, increasing the aperture size, the width 
of the Gaussian function decreases, reducing the dip vis- 
ibility (see Fig. [2]). Physically, this can be explained 
by the fact that increasing the aperture size we let more 
wave-vectors into the system, and so the spatial walk- 
off in type-II interferometry introduces distinguishability, 
reducing the interference visibility. Enlarging the aper- 
ture sizes is often useful in practice, for example to get a 
higher photon flux. Moreover, since in the SPDC process 
different frequency bands are emitted at different angles 
it may be necessary to open the detection aperture in ap- 
plications where a broader bandwidth is needed. This is 
clearly a problem when using type-II phase matching in 
birefringent crystals, since the visibility of temporal and 
polarization interference gets drastically reduced. 



To get an analytic result one may assume Gaussian 
detection apertures of radius Rq centered along the sys- 
tem's optical axis: 



B. Linear phase shift 

Suppose now we introduce a linear phase function with 
the spatial light modulator, along the directon Si: 



p(x) 



(29) 



Si • X 



(34) 
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we get: 



Wm{t) = Sine 



Pa 





( 2r M 




2diD 





Mkp 
2diD 



Te-2 + 2/si 



Pn 



Mk 



■p 

Tea 



2/si 
(35) 



Therefore, the amount of tilt necessary to restore high- 
visibiUty is: 



tan 6 — 



ML 



(40) 



In the case of a 1.5 mm crystal, with M = 0.0723 (pump 
at 405 nm, SPDC at 810 nm) and lenses of focal length 
20 cm in the 4-f system, we get: 



0.14 mrad 



(41) 




T(fs) 



FIG. 3: Quantum interference pattern for different detec- 
tor aperture sizes, introducing a linear modulation of the de- 
formable mirror, in order to restore the indistinguishability 
between the photons, decreased by the spatial walkoff in the 
generation process. 



If we compare Eq. (|35|) with Eq. ([30)1 we can see that 
the structure is the same; we again have a triangular 
function centered in r = DL/2, along with two aperture 
functions. But this time, instead of being centered at 
T = 0, the aperture functions can be shifted at will along 
the T axis. Suppose we now apply a tilt along the y axis 
{■Six = 0); the modulation function is then shifted to: 



^cente 



ID 
koM 



Sly 



(36) 



To get the highest possible visibility, the center of the 
modulation function must be matched to the center of 
the triangular dip: 



so that: 



^center 



DL 
k^ML 



(37) 



(38) 



In the case of a reflective system, in which the phase 
modulation is implemented by means of a deformable 
mirror, tilted by an angle 9: 



C. Large aperture approximation 

If the detection apertures are large enough for the 
Pi function to be successfully approximated by a delta- 
function we get: 



(42) 



Suppose that the spatial modulator is a circular aper- 
ture with radius r, with unit transmission and phase mod- 
ulation described by the function V'(x): 



r if |x| > r 



(43) 



In this case the function <p(x) can be expanded on a set 
of polynomials which are orthogonal on the unit circle, 
like the Zernicke polynomials: 



<^(q) = 51 XI "P^mRnip) C08{me) 
n m 

(m = — n, —n + 2, —n + 4, n) 



(44) 



(/5(x) = 2fco tan^y = siy-y 



(39) 



where q = {pcos9, psm9). To calculate ip{—q} we note 
that — q = [pcos(6' + n),p am{0 + tt)], so: 

^i-^)^T.T.^nip)^os[m{0 + 7r)] (45) 



If m is even then coa[m(9 + tt)] = cos{m9), otherwise if 
m is odd cos[m(9 + n)] — ~ cosni9). Therefore: 

<^(q)-V'(-q) E ^nmK{p)cos{m9) (46) 
" m odd 

So, only the Zernike polynomials with m odd con- 
tribute to the shape of the interference pattern. This 
effect is the spatial counterpart of the dispersion cancel- 
lation effect, in which only the odd-order terms in the 
Taylor expansion of the spectral phase survive. The ex- 
perimental demonstration of this effect has been recently 
reported in |11|] 
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FIG. 4: Example of the numerical approach adopted to evalu- 
ate Eq. (|25|l and Eq. (I26p . The spatial modulation surface is 
discretized in sufficiently small squares over which the phase 
is averaged. 



IV. NUMERICAL SOLUTIONS 

Numerically solving for the quantities in Eq. (|25p and 
Eq. (j26p may be computationally demanding. Here, we 
propose an approximation, valid in the case where the 
function G(x) changes smoothly over the mirror surface, 
as it is the case in experimentally relevant situations. 
This model is also interesting from the practical point 
of view, since in many cases adaptive optical systems are 
implemented using spatial light modulators or segmented 
deformable mirrors, where the modulation surface is par- 



titioned into small pixels. 

Suppose we partition the Fourier plane Hp into small 
squares (pixels) of side d. Let's define the rectangular 
function: 



n(x) 



if |x| > 

1 if Ixl < 



(47) 



The pixel (1, m) is identified by: 

X 



'^i,mix,y) = n 
selecting the area: 



n 



(48) 



(49) 



We approximate the value of the phase in each square 
by the mean value of the actual phase within the square: 



<Pim = I dxdyip{x, y)II 



n 



y 

- + m 
a 



that is: 



(50) 



(51) 



In this case (see Appendix B for a justification): 



+ 1 



n 



Lm 



l,rn 



(52) 

Substituting this expression in Eq. 1231 and collecting the 
integrations one finds: 



l,m X,fi 



'■'^aixl^^ir) (53) 



where: 



dqx / dQxU 



^ 1 



n 



(54) 



and: 



Im,ti{T) = / dqy / dQyIl 



f_ 



n 



Sine 



ML{qy + Qy)A[l-^ 



(55) 



P[qy + Qy]P[-{qy + Qy)] 
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Performing the integrations one gets: 



dxP{x)K 



fx 
kd 



(l + X) 



2ddi . 
sine ■< — —xA 



fx 
kd 



(l + X) 



(56) 



and 



W(r) / dxP{x)A 



fx , s 



sine < MLxK 



{2kdf2di M 
inr, < X 



I / 



D 



T A 



1 



fx , ^ 



2r ' 

Ul 



(57) 



A similar expression can be found for the background 
coincidence rate: 



l,m X.fi 



(58) 



where: 



(x) 
IX 



dxP{x)A 
2ddi 



smc 



f_x 
kd 



xK 



f_x 

kd 



-il-X) 



e f 



(59) 



and 



dxP{x)A 
2ddi 



fx 
kd 

xA 



[m — fi) 



sine (MLx) 



fx 
kd 



(60) 



The advantage of our numerical approach is that one 

fa:) (v) 

can calculate and tabulate the functions , i?m/i, ctix 
and Imn{T) for a given configuration, determined by the 
focal length f, the shape of the detection apertures, the 
width of the deformable optics and the distance between 
the crystal and the detectors. Then, to calculate the 
shape of the interference pattern for a certain phase dis- 
tribution on the adaptive optics one just needs to change 
the coefficient of a linear combination of the tabulated 
functions. This can be a helpful tool to study the effect 
of specific aberrations on the temporal interference or to 
engineer the shape of the HOM dip. 

In the limit of large detector apertures P{x) « 5{x): 



a/A = 5{l + \) 



Imn = sine 



2kdM 



(61) 



(62) 



So: 



Rbig{t) « E' 



-iiVlm-iP-l.-m) 



ai,_lIra-m{T) (63) 



We have a discrete formulation of the even-order spatial 
phase cancellation effect. What affects the shape of the 
dip is the difference between the phase of the pixel (1, m) 
and the phase of the pixel symmetric to it with respect 
to the origin (-1, -m). Therefore, since for even func- 
tions (pim = (j)-i^-rm the phasc difference is zero and they 
do not contribute to the coincidentce pattern. This way 
our approximate solution technique is consistent with the 
general theory. 



V. DISCUSSION 
A. Examples 

According to Eq. [531 the shape of the temporal quan- 
tum interference pattern is affected by the spatial mod- 
ulation given by the deformable mirror. In Section IV 
we proposed a numerical approach to calculate easily the 
shape of the temporal interferogram. For example, some 
results are reported in Fig. [5l for coma (upper plot) and a 
superposition of several different aberrations (lower plot). 
The interference visibility clearly degrades in presence of 
wave-front aberrations. 



B. Limits and role of the large aperture 
approximation 

An interesting question is under what experimental 
conditions the large aperture approximation can be con- 
sidered to be valid, so as to obtain the even-order aber- 
ration cancellation effect. 

According to the numerical approach proposed in Sec- 
tion IV, the even-order aberration cancellation effect 
manifests itself in the limit where, for example: 



ai,x S{l + X) 



(64) 



In Fig. [HI a plot of the value for ai^\ is shown for 
typical values of the relevant experimental parameters 
(detection aperture radius R = 5 mm, detection distance 
di = 1 m and size d = 0.1 mm of each pixel in the Fourier 
plane of the adaptive optical system). Clearly, only the 
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FIG. 6: Plot of the values of |a!,Ap, for /,A = -100... + 100. 
The radius of the detection apertures is _R = 5 mm, the dis- 
tance between the exit plane of the modulation section and 
the detection apertures is d — 1 = 1 m and the size of the mod- 
ulation pixels is d = 25/im. Clearly only the diagonal elements 
contribute are non-zero, i.e. the ones for which A — ~l. In 
this situation the effect of even-order aberration cancellation 
is present. 



200 400 

tCfs) 



FIG. 5: Examples of the shape of the polarization quantum- 
interference dip with two different spatial phase modulations 
(in black the unperturbed dip, in red the modulated one). In 
the upper figure the effect of a small amount of coma along the 
vertical axis is shown. In the lower figure a more complicated 
superposition of aberrations affects dramatically the shape of 
the dip 



diagonal elements (the ones for which / = — A are signifi- 
cant), suggesting that the effect of even-order aberration 
cancellation may be observable for most typical experi- 
mental parameters. 

To get an idea of what happens for different experi- 
mental conditions we can compute the ratio between the 
intensities of the non-diagonal coefficient aoi and the di- 
agonal coefficient aoo: 



Foil 
I aoo I 



(65) 



The lower the value for pQ the less significant the coeffi- 
cients for \ ^ —I are: the even-order aberration cancel- 
lation effect will therefore manifest itself more clearly. 

Values for po are shown in Fig. [7] for two different cases. 
In both pictures, the value of po is shown as a function 
of the Gaussian detection aperture radius R, for three 
different values of the distance between the plane li^ and 
the detection lenses di. In the upper figure, the size of 
each small square in which the spatial phase is assumed 
to be constant is d = 0.5 mm, while in the second case it 



is d = 0.25 mm. In both cases po is significantly smaller 
than 1, and it becomes smaller and smaller increasing 
the value of the detection aperture radius. However, po 
is smaller for larger values of d, implying that the spatial 
variability of the modulation phase plays a role in the 
degree of even-order cancellation of the modulation itself. 



It turns out that for the aberration cancellation effect 
to appear, it is in fact only necessary for one aperture 
to be large and for one detector to be integrated over. 
This is sufficient to produce the transverse-momentum 
delta functions that lead to even-order cancellation. To 
demonstrate this, we can for example consider the case 
where the aperture at B is large, and the detector at B 
is integrated over, while the aperture at A is taken to 
be finite, with detector A treated as pointlike. The lo- 
cation of the pointlike detector will henceforth simply be 
denoted as cc, and we continue to work within the quasi- 
monochromatic approximation. If we integrate only over 
xb, leaving x unintegrated, then it is straightforward to 
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Fourier transform of |p_Bp. We now let aperture B be- 
come large, so that the function Pg goes over to a delta 
function. For Gi{x) = e'^'f'(^) and G2{x) = 1, we can 
substitute these results into the coincidence rate (which 
will now be a function of both r and the position x of 
detector A), and carry out the q' and v integrals. For the 
modulation term, wc find: 



Rm{x,t) = R{x,t) ~ Ro{x) 

= I dq e^[^(9)-^(-?)l 



(68) 
(69) 



2i,2 r 


4f+i 






c/o y 


cfo 





Sine 



c/o 



Here, we have used the fact that the Fourier transform of 
p\{x) equals the complex conjugate of the Fourier trans- 
form of PA (a;), in order to write Q' in terms of Q. We 
see from the presence of the factor e*^'^^''^"'^'^"*^! that 
even-order aberration cancellation occurs even though 
one aperture is finite and the corresponding detector is 
pointlike. This point may be of importance in future 
attempts to produce aberration free imaging. 



FIG. 7: Plot of the ratio po between the intensities of the 
non-diagonal coefficient aoi and the diagonal coefficient aoo 
as a function of the radius of Gaussian detection apertures 
R, for different values of the distance between the exit plane 
of the modulation section and the detection apertures (di = 
1, 10, 100 cm). On the upper plot the size of the modulation 
pixel d is d = 0.5 mm, while in the second case it is d = 0.25 
mm. Clearly, for experimentally interesting cases, the off- 
diagonal coefficient is at least 3 orders of magnitude smaller 
than the diagonal one, leading to even-order aberration can- 
cellation. 



show that: 



(66) 



(67) 



-q + — r ]Q +9 + —r ]PB{q-q) 

cfo J \ cfo J 

Here, we have defined Q and Q' to be the Fourier trans- 
forms respectively of pA and p*^. Pb is, as before, the 



VI. CONCLUSIONS 

Summarizing, we have in this paper done a theoretical 
study of the relation between the wavelength modulation 
of the entangled SPDC photons and the shape of the re- 
sulting temporal quantum interference pattern. Due to 
the multiparametric nature of the generated entangled 
states, the modulation on the spatial degree of freedom 
can affect the shape of the polarization-temporal interfer- 
ence pattern in the coincidence rate. Our aim is twofold: 
from one side we want to study the effect of wavefront 
aberration on quantum interferometry, from the other we 
want to discuss a way to engineer multiparametrically- 
entangled states. 

We have introduced a theoretical model to calculate 
the shape of the polarization-temporal interference pat- 
tern given a certain phase modulation in the crystal far- 
field, assuming as a free parameter the shape and the 
dimension of the collection apertures. Using a numer- 
ical method to study the resulting equation has shown 
that for typical experimental cases the hypothesis of large 
apertures can be assumed valid. In such an approxima- 
tion, only the odd part of the assigned phase modulation 
affects the shape of the interference pattern. This effect 
has recently been demonstrated experimentally [111] . 

Moreover, it is often useful in experiments to enlarge 
the collection aperture in order to collect a higher photon 
fiux and larger optical bandwidth. But when working 
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with type-II birefringently phase-matched downconver- ([16]), and the result into Eq. (|22p one finds the following 
sion, spatial walk-off between the emitted photons intro- expressions for Rq and WciT)'. 
duces distinguishability between the two possible events 
that can lead to coincidence detection, reducing the visi- 
bility of quantum interference. Such walk-off can be com- 
pensated with a linear phase shift in the vertical direc- 
tion, restoring high visibility. 



XG* Uq\ G2 M^W(q,q',z.) 
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APPENDIX A: DERIVATION OF EQ. 



(A2) 



In this Appendix we sketch the major steps for the 
derivation of Eq. (^5)) . Substituting Eq. (PT|) into Eq. where: 



and 



M^(°)(q,q',i/) - J dxAdxBi/*(xA,q, J^)if*(xB, -q, -i/)Jy(xA, q', (xb, -q', -i^) 

+ -ff*(xA, -q, -i^)-ff*(xB, q, i^)h{xa, -q', -i^)h{xb,c[', v) 

W{(l,c{,u) = y dxAdxB-ff*(xA,q, i^)iJ*(xB,-q, -J^)i?(xA,-q',J^)-ff(xB,q',-!^) 
+ H*{xA, -q, -J^)-ff*(xB, q, J^)-ff(xA, q', -i')-ff(xB, -q', i^) 

I 



(A3) 



(A4) 



The angular and spectral emission function $(q, j/) is and 
given by: 



1 dz n 


' z 1" 






L + 2 



-■iA(q,f)2: 



(A5) 



Performing the integrals over the spatial coordinates dxA 
and xb one gets: 



I^(q, q', u) = e'^[l''l'-l'''l'lpA[(q + q')]^^i3h(q + q') 



PA[-(q + q')]^i3[(q + q')] 



(A7) 



W^('')(q,q',j.) =e^T^[l''l -I'l'I'lp^Kq - q')]PB[-(q - q')] 
+ P^[-(q-q')]Pi3[(q-q')] 



Finally, use of the integral representation for the sine 
function (Eq. (|A5p ) allows the v integration to be carried 
out, but at the expense of introducing two integrations 
(A6) over a pair of new parameters (say z and z'). Note the 
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relation 



If we express the first few terms we get: 



1 ,if-l<a<0, -i<a;<i + a 



n[a;]n[x-a] ^ { 1 ,ifO<a<l, -i+a<a;<i 
, else . 



From this, it follows that 



j n[a;]n[a; - a]dx = A(a), 



(A8) 



(A9) 



where A{a) is the triangle function. These facts allow us 
to carry out the two z-integrations that arise from the 
sine function, leading to the result shown in Eq. (|23p. 



APPENDIX B: JUSTIFICATION OF EQ. ((521) 

Suppose to have a set A, which can be partitioned into 
a collection of disjoint subsets Ak, k = 1,2, 



[jAk^A AkCiAi^c^ if k^l 



(Bl) 



fe 

= [1 + (e**^ - l)xi] [l + (e'^^-l)x2]... 

= l + (e'^i -l)Xl + (e^^^-l)X2 + ... 
+ {e'*' - l)(e^*^ - 1)X1X2 + 

+ (e'^^^-l)(e^*^-l)xiX3 + .... 

+ (e'0i _ i)(e^0i - l)(e»^i _ l)xiX2X3 
+ (e*0i - l)(e^<^^ - l)(e*^^ - l)xiX2X4 + 



So that in the end: 



(B6) 



To each set we can associate a characteristic function: 



Xk{x) = 



1, X e Ak 



(B2) 



_ 0, x^Ak 
such that: 



where xa is the characteristic function for the full set, 
XAix)^{l^ Hi (B4) 



The term e^'^''^'^'^^^ assumes the value e"^'= for Xkix) = 1 
and the value 1 for Xfc(a:) = (1 — Xk{x) = 1), so: 



Y[[l-il-Xk{x))+e^'^''-Xkix)] (B5) 



giE. 'PkXkix) = 1 + ^ [(e'-^-^ - 1) xfe] 

k 

= l + Y.''''Xk-J2xk (B7) 



Xk 



Since the square sets we have used in section IV satisfy 
Eq. (Bl), then the result expressed in Eq. (B5) is valid 
for our case. 
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